' 1 - ' y ' j 
■'rxiAOO 

• y> a .'A' ~ V 
> X ; \TVl/ 


yyA'AX 

■ V • "* 

5 • \ 

GPO PRICE $ 

/ - 

CSFTl PRICE(S) $ 

V 


X-'X' , \/ , 

4, ; ' , 

XAj&AyK v 'X ■ 
, ,R xy ; ;. H 


' : ; v ;. ;v v !i? '/ >■', - z ' A'; .’ax ■: .\X - a a/> : :a’' ax»a : ;< a 

.' ': i-,'-' v;r' .7v-,\ '"v. • A' / v A' v , ■ " ■ 

* *•• - - i’(S iX ;• ! k) K; />y VX S'Y'z : V'/;)/- ; X-643-fi7-6t3 V 

■ Vi- : \ ^ , Vu ,d xf ' V £ \V_*\ • x / V p .iy yl\ X.rx vx //j V -,X 

' f '‘‘n'V ■ ‘ f ' v V.'ft v ,\ ) !'vi: * - 4 l ) / ■’. l A ; ••.;• (■ . '-ffyvt .,1" I' 'X. 

:x v '^ j'x J/x 1 / S X‘ : n>i' (;/' \>\x *U\i ' ,'y *x , .• H V V .^; v- 

\> : r . 

► Y ’w' v. : A ’ «■ / t ' ' . ,{ / ii f f \ x „ •! v~ , - | lx i A , x< v/; • *j , 

% $ W)M NEW POWER SERIESINVERSIONS 

XyX i'V" .■ ' ' ; '-'p; x v. «!.' ';', ■ .'. ',y lv :s; -^ v«, 

,.f : ■ v r ;\h. ;•■ 

• v'K s I ,; !'h; 1 < v v /Y '-A , ; ^' c, xY\ , 

V( ,} , , ^ ( * 1 J ** ^ ^ / Y' ? %4';’> ••''•/% •* •, .V\ ' /, i ■ 

"-.V .': ( 'Y A; Vu-Y ' ': AT ' ' ■■ "' ; f • v l; v' r ’ i ! , ^ '4 : |X-: A -i r - ’ 

' ' </'■ v X " : s '<v a:' ' >/■; V?' •''?•■■ - 'r A ■ Q 

* h/ ' ' ■ ■ h • !y.:f! ^ ■$ ■ Eu' ^ 

40 uO\$: rxA'xO Hi 

0004^04 

nOAOOOf 

J-J: '■ . v ) A O- A 'i i ■’ 4 x ' [ x 

; aY 5 vfilX 'A ;X :r x -> 


Microfiche (MF) _ 


j ff 653 July 65 

(\'v ■ •■xxy- % ) 

AfOO 

A y4- Ax. 

■ i "4 AA 'O ■ 

4‘k0A • 

'■>o- }o\h,. 

xx A .A' '■-■?<: 

’T \ v Kufy v 

• a \ -*4 

*i 040. 0 

v ' '>y„ -v 


AOaooAoa 

lAOOAO.p 

%Oo : o 

-x;¥\ ):n',f 

A-v'v" 


k J 


;trf 


‘^•V- ’■ \ )r K X- i 

... A\4 \ 


■M A , 



(ACCESSION NUMBER) 

12. 


(THRU) 


(PAGES) 


2 

£ 

_I 

^ (NA f SA CR OR TMS OR AD N^mSeR) 


(CODE^ 


z 




(CATEGORY) 


> vv x,v • • :,. • y ^ , y ■ A’. ( ■ 

<4^. -a ’ • ’ - **-sJ 


YU ' ’V " ■ >fp; hkI\- : i' >•* 'f" ; 

, \ > M ’'' y '/X -•'/ A ' A V'v 1 v .j 

V * », ./ - Th\ \ -'\v .*&..« - -. , ^fXllv7 T '''kXW i X" ';:'V.i : 'X>. (/X Aa^a,vV;VX- /yiA^: t y z Mz 

• V ; . A ; '^ j 'i’.CjVi ' ; • J; 1 'l >>,\A • -YV'V''' 9 ' A ^ v '/ ' 

J- ■; ,' ;M - Aik ixXxiX’Xx 


■i, t :r- t\ r > ■ v r - x^ r 40 WApiA " AAA ; .< y ; , 



ax -A, yr -5 z 4 . ro- 1 ' 1 


o n ^ 

ffAXA ) 

AA A; - x y ' \ V ^ <i (:. 


' V J A 3t\V: j ■ ’W 


‘ 1 •*;. •’. ■ • ;; i ;x'‘ x i>u x i >. a 'Qa':' ‘ c , 

, ' x , > i • ' V v .V 'A : A \ 4 A 'A ' ' V.’ 'j.Y'.j' . .A \j! i; ^' A ' ; y', ',!/ ! i \ '.'J--' 

-X; #kf*a«%F ; : "w lAiiT utfii 


A A .A ’ M v / ' ' ! 


1 i'V 



f - 


'4XA 


^ . X 


^Wm^A: 'v MARYLAND 


'' " ' Ai ■-' . ‘ 


■/. ' . - K < 


■•■ v^: X.Aj V x :,i 

\i AO ■>.-. 4 'if. • '• • X.\r ': A. 0 ■ y 


?| , ^ ' 




\y,- v;i - 

•. t -v 


M 


■j > 


> X/ 


1/ .. ! , •.•■• .-A V X- . .. . s%, V ; . ■ • ', 

■ 'y . ••. y \ x , \ y ' , > j :•. \, . y 

- -x^ -y'\ V : ■* ,/y '*' ,; y y 


Sr; -P - a; i >; \ o - X -- ]{ < x :y 

fAAOfOo 4 AA 04 AX. 44 AA 

I XfXl4 . Vi i'OA-y : '.''•l' ;" ^y 

■• ; XVy '1 A OA ;•' > j \ 0 r^^ ' v A : ' 

'. II i " zz wOfi 1 ! A ^M/OOAAi'H. X 

AO^A?i4lOOOi4A4 

, X' ! •./— oj4.yJ yyxAy.' 

i’i\ . ! 'y ?9 ’-:f*y /i ’ -1 - r > :A 



-1.,: ,J_ J_ ,:AA ii .::I AAA...:. .A. 


X-643-67-613 


SOME NEW POWER SERIES INVERSIONS 
OF KEPLER'S EQUATION 


E. R. Lancaster 
R, H. Estes 


December 1967 


Goddard Space Flight Center 
Greenbelt, Maryland 




LIST OF SYMBOLS 


a semi- major axis of orbit 

e eccentricity of orbit 

E eccentric anomaly of orbit 

E 0 eccentric anomaly at time t 

M mean anomaly 

r magnitude of the orbital position vector 

r 0 magnitude of the orbital position vector at time 

r time derivative of r 

r Q time derivative of r at time t 

t time 

/x gravitational constant times the mass of the attracting body 

x component of the position vector on an axis with origin at the 
primary body and positive direction toward pericenter 

y component of the position vector on an axis in the plane of 
motion 90° from the x-axis in the direction of motion 





ABSTRACT 

Recursion formulas are derived for the coef- 
ficients of three power series inverses of Kepler's 
equation: 

1. a series in powers of time, 

2. a series in powers of eccentricity, 

3. a bivariate series in powers of two small pa- 
rameters. 


For each series a formula is given for the radius 
of convergence. 



SOME NEW POWER SERIES INVERSIONS 
OF KEPLER’S EQUATION 

THE CLASSICAL SERIES 

Let E be defined as a function of e by Kepler’s equation 


E-esinE-M. (1) 

We wish to express E as a Maclaurin series. This is customarily done [1] 
by means of a Lagrange expansion. The procedure we follow is simpler and 
produces recursion formulas for the coefficients in the series. Moreover, the 
method simultaneously produces series for sin E and cos E. 

Let 


u = sin E, 


w = cos E. 


Then 


u* =wE\ 

(2) 

w ' = -uE\ 

(3) 


where the prime indicates differentiation with respect to e. From (1) we also 
have 


E - e u = M. 


Let 


E = 



a. 

1 


(4) 


(5) 


1 



Then 


uu 

u=2] b i ei > 


£c. ei - 


1 ‘ = Tl (i + 1} a ' i+i ei> 


uu 

“• = L <* 


i + 1) b i+1 e\ 


(i + 1) c. e\ 


Substitute (5) through (10) in (2), (3), and (4) and equate coefficients of like 
powers of e to obtain 


a o =M > 


b Q = sin M, 


c Q = cos M, 


a. +1 =b., i *o, 


d , =3 


2 



(i + l)b i+1 = d j+ i c i_j' 1 “°> 

J^o 


(i +1)Ci+i= Z] d j+ ibi -i’ i “°' 

j=0 

Having series for sin E and cos E, we can easily obtain series for r, x, 
and y from the following formulas. 

r = a(l - e cos E), 

x ~ a (cos E - e), 

y = (1 - e 2 ) 1/2 s in E . 

Plummer [2] uses Lagrange's theorem to obtain the radius of convergence 
for these series. We have obtained the same result by the following simpler 
procedure. 

The radius of convergence is equal to the smallest magnitude of the values 
of e associated with the singularities of the function E (e,M ), which occur where 
E' ( e, M) is infinite. Since 

v, _ sin E 

- — r 

1 - e cos E 


the singularities are determined by 

1 - e cos E = 0. 

Substituting this in (1), we have 

E - tan E = M. (11) 
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Let E = u + iv and E = u - iv. Then 


tan E = 


sin E cos E _ sin 2 u + i sinh 2 v 
sin E cos E cos 2 u + cosh 2 v 


( 12 ) 


Let w =t I e | 2 , Then 


w =T e e = 


cos 2 u + cosh 2 u 


(13) 


Substituting (12) in (11) and equating real and imaginary parts gives 


u - w s in 2 u - M , 

(14) 

v - w sinh 2v = 0. 

(15) 


Equations (13), (14), and (15) give u , v, and w as functions of M. We de- 
sire the values of M which minimize and maximize | e | or, equivalently, w. We 
are thus led to set dw/dM = 0, obtaining 


~ sinh 2v - sin 2u = 0 , 
dM dM 


— - (1 - 2 w cosh 2v) = 0, 


du 

dM 


(1 - 2 w cos 2u) = 


1 . 


These equations imply sin 2u = 0 or cos 2u = ±1. Substituting this in (15), we 
obtain 


v cosh 2v = sinh 2v ± v, 
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which can be reduced to 


coth v = v 


(16) 


for the + sign and to 

tanh v = v (17) 

for the - sign. The solution of (17) is v = 0. The solution of (16) is approxi- 
mately 1.20. Substituting these values in (13) we obtain | e ] - 1 and ] e | = 0.663 
for the maximum and minimum radii of convergence. 


A BIVARIATE SERIES 

For M in equation (1) we have 

M = (t “ t 0 ) (Av'a 3 ) 1/2 + E 0 - e sin E Q . 

Because E 0 in this expression is indeterminate for circular orbits, Kepler's 
equation is often written in the form [3] 

y — a sin y + fi (1 - cos y) = r, (18) 

where 


7 = E-E 0 , 


a = e cos E Q = 1 - 


P = e sin E 0 = r Q r 0 /(/^a) 


1/2 


T = (t - t 0 ) (/Va 3 ) 


3 \ 1/2 


Since a and fi are small when e is small, it seems natural to seek a Maclaurin 
series in a and fi. 


5 



Let 


u = sin y, 


w = cos y ■ 


Then 


7 - au - /3w + /3 = T , 

(19) 

3u 3 7 

3a ~ W Ba 

(20) 

3u 37 

3/3 ~ W 3/3’ 

(21) 

3w _ 37 

3 a 3 a 

(22) 

3 w 37 

3/3 U 3~/3* 

(23) 

00 °° 

r= Z] Z] an* 1 ' 31 - 

issO j srO 

(24) 

00 00 

•• LL 

issO jssO 

(25) 

00 00 

w= Z! Z] 

issO j ssO 

(26) 


Substituting (24) through (26) into (19) through (23), using the formula derived 
in the appendix, and equating coefficients of like powers, we obtain the following 
recursion formulas. 
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a oo = T ’ b oo = sin r, c 00 = cos r 


a io ~ b oo ’ a oi ~ e oo * 


i + 1 ,0 ~ 

c r 
►**• 
o 

a o , j +1 =c 

3 i+l . j 

= b i3 

+ C i+1, j-l 

a i.3+l 

= b u 

1. j+1 + C ij 

PC 

II 

(m + 

a m+l , k 

& 

3 

PC 

n 

(k + 

a m,k+l 


Oj 


3 i 


(i + 1) b 


i+l,j ^ ^mk C i-m,j-.k 

k— 0 m-0 

3 i 

0 + 1 ) b i, j+ l= 2] £ qmk C i-m,j-k 

ks 0 niss 0 


J 1 


(i + !> C i + l,j J2 H P mk b i-m, j-k 


kss 0 m=0 


J i 


( j + 1) c i ^ j + i * lz: ^mk b i-m , j -k 


k-sO m=0 


These series converge for 


a 2 + p 2 < 2 w , 

where w is found as for the classical series. When computing with these series 
the coefficients are found in sets such that for the nth set the sum of the sub- 
scripts on each coefficient is n. 
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A TIME SERIES INVERSION 


If we consider y as a function of r in equation (18), we 
Maclaurin series inverse. Letting u = sin 7 and w = cos 7 

can again obtain a 
as before, we have 

7 - au - /Sw + /3 = T , 

(27) 

u = w 7 , 

(28) 

w = - u 7 , 

(29) 

where the dot indicates differentiation with respect to r . 


Let 


CO 

7 - a. r 1 , 

i = 0 

(30) 

00 

- 

i = 0 

(31) 

00 

w = ^ * c i r 1 . 

i=0 

(32) 

Substituting equations (30) through (32) in equations (27) through (29) and 
equating coefficients of like powers of r yields 

c 0 - 1 , b Q - a Q - c x 0 , 


d = a/r Q , = d 


C 2 ~ ~ ~2 ^.1 ’ k.2 = .fi ^ C 2 ’ a 2 = ^2 


= J a. 
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C i+1 _ " 


b i + l = d 


22 7 > b i-j t i + i) 

^ c i + i + (^22 7 > ^U-» 


a . , = a b . - + >S c . t 

i+i l+l ^ i+i 


It is not difficult to show that the radius of convergence for these series is 


r - 



ar ctan 



a 


)* 


+ (sech -1 e - /l - e 2 ) 2 


1/2 


where 


-7r < a r ctan — 5 77 -. 

a 


COMPUTATIONAL CONSIDERATIONS 

If one desires merely to compute updated values of E or y and associated 
quantities, the following alterations in the algorithms reduce the amount of 
computation required. 

For the classical series replace the equation a i+ i = b i > 0 by 

a i+ i = e b i • i - 0 

so that the resulting series then become 


E = 



iatO 


9 



s in E - 


CO 



i=0 


For the bivariate series use 

a i 0 = a b 00 ’ a 01 = ft ( C 0G “ 

a i + i,0 =ab i0> a 0 , j +1 =^ C 0J 
a Hl,i = ab ij + / 3c i + l,j-l 
a i,j + l =ab UX,j + l + ^ c ij 
The series then become 



Finally, for the time series we make only the change 

bj = a i = r d, 
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and the series become 


7 = 



sin 7 = 



cos 7 


00 

" 2 > 

isO 
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APPENDIX 


BIVARIATE CAUCHY PRODUCTS 


The Cauchy product of the series 


S 


1 


E 


a. z 


and 


S 2 = 


CO 

E 

1=0 


b. z 


is defined as the series 


S = 



c i z 


where 


i 



j 5=0 


We extend this definition in a natural way to the product of two bivariate series. 
Let 


00 00 


EE a ‘ ial/3i = Z> Bi - A i = E a “ 

j = 0 i=0 j = 0 i=0 


12 



s 


2 


00 00 



}=0 i = 0 


b.. a 1 fii = ^ B j B i = 

j «0 


00 



iasO 


Applying the above definition of a Cauchy product to these two series 
gives 


s = s .Sr £ 

j a0 



CD 00 



jatO i = 0 


c ii 



kssO m=0 
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